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ABSTRACT 

The flux-ratio anomahes observed in multiply-lensed quasar images are most plau- 
sibly explained as the result of perturbing structures superposed on the underlying 
smooth matter distribution of the primary lens. The cold dark matter cosmological 
model predicts that a large number of substructures should survive inside larger halos 
but, surprisingly, this population alone has been shown to be insufficient to explain 
the observed distribution of the flux ratios of quasar's multiple images. Other halos 
(and their subhalos) projected along the line of sight to the primary lens have been 
considered as additional sources of perturbation. In this work, we use ray tracing 
through the Millennium II simulation to investigate the importance of projection ef- 
fects due to halos and subhalos of mass m > lO^/i^^M© and extend our analysis to 
lower masses, m ^ lO^/i^^M©, using Monte-Carlo halo distributions. We find that the 
magnitude of the violation depends strongly on the density profile and concentration 
of the intervening halos, but clustering plays only a minor role. For a typical lensing 
geometry (lens at redshift 0.6 and source at redshift 2), background haloes (behind 
the main lens) are more likely to cause a violation than foreground halos. We conclude 
that line-of-sight structures can be as important as intrinsic substructures in causing 
flux-ratio anomalies. The combined effect of perturbing structures within the lens and 
along the line of sight in the ACDM universe results in a cusp-violation probability 
of 20-30%. This alleviates the discrepancy between models and current data, but a 
larger observational sample is required for a stronger test of the theory. 

Key words: gravitational lensing: strong - galaxies: haloes - galaxies: structure - 
cosmology: theory - dark matter. 



1 INTRODUCTION 

In the cold dark matter (CDM) cosmogony, galaxies are bi- 
ased tracers of a filamentary "cosmic web" of collapsed re- 
gions in the matter density field - dark matter haloes. The 
excellent agreement between the predictions of this model 
and observations of the large-scale clustering of galaxies pro- 
vides compelling support for CDM. However, on the scale 
of individual dark haloes, the model makes a number of 
predictions that have yet to be fully verified: cuspy halo 
density profiles and a large population of surviving dark 
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matter substructures. These substructures are the cores of 
accreted CDM haloes that persist as lon^-lived gravitation- 
ally bound subhaloe s (jCao et al.ll2004l : lDiemand et al.ll2008l : 
'Springel e t alll2008l ). Therefore, measurements of halo den- 
sity profiles and of the subhalo abundance are crucial tests 
of the cosmological model. 

Galaxies and their dark matter haloes can act as strong 
gravitational lenses, producing distorted and even multiple 
images of more distant galaxies and quasars. The distribu- 
tion and properties of these images provide sensitive probes 
of the mass distribution in the lens. In some multiply-lensed 
quasar systems, simple parametric mass models can fit im- 
age positions well, but not their flux ratios. Such anoma- 
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lies are interpreted as evidence for complex s ubstructures in 
lensing galaxie£] (Mao & Schneider 1998; M etcalf k Madaul 
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2001; Metcalf & Zhao 2002; Chiba 2002; Me tcalf et al 
Sugai et al. 20 07; McKean et al. 2007; More et all [200! 
MacLeod et"al]l2009l ). 

On scales probed by galactic strong lensing (typ- 
ically a few kiloparsecs) , predictions from CDM sim- 
ulations have been compar ed with observed flux-ratio 
anomalies fe.p:. | Dalal Kochanekll2002l : iBradac et al. 2004 : 
iMetcalf k Amaral l2012l C Several studies have concluded 
that the predicted abundance of dark matter substructures 
in the strong-lensing regions of galaxy-sized haloes is not 
sufficient to explain the statistics of the currently avail- 
able sample of known anomalous lense s (|Mao et al. l2004l: 
Amara et all l2006l : iMaccio Mirandal l2006l : IChen et al.1 
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In any smooth lens potential producing multiple im- 
ages of a single source, a specific magnification ratio (here 
equivalent to a flux ratio) of the three most strongly- 
magnified images will approach zero asymptotically as the 
source approaches a cusp of the tangential caustic. This 
is known as t he "cusp-caustic r elation" (see Eg . p be - 
low) (iBlandford k Naravanlll986l: ISchneider k Weissl]l992 



IZakharovl Il995l : iKeeton et al.1 \20oK ~ Structures, either 
within a lensing galaxy or projected by chance along the 
line of sight, will perturb the potential and alter the flux 
of one or more images, resulting in a violation of the cusp- 
caustic relation. These violations are extreme cases of flux- 
rati o anomali es. 

Xu et al. 1 (|2009l . I2OIOI ) analyzed flux-ratio distributions 
of multiple-imaged background quasars in simulated lens- 
ing systems, using six ^ 10^^ Mq CDM haloes and their 
substructure pop ulations (subhaloes a nd streams) from the 
Aquarius project ("Springel et aLl fioOSI V The effects of bary- 
onic substructures (satellite galaxies and globular clusters) 
were also investigated. These exceptionally high resolution 
simulations confirmed that the substructure abundance (in 
the critical region of a Milky Way-mass lens) predicted by 
the CDM model is too low to explain the observed frequency 
of cusp-caustic violations. 

This apparent deficiency of substructures is not yet 
a strong challenge to CDM, in part because the sample 
of observed lenses is extremely small. Furthermore, dark 
matter haloes and subhaloes are present along the entire 
line of sight from the source to the observer, not just in 
the lens itself. These independent haloes projected in front 
of and behind the lens can also induce perturbations to 
the lensing potential and thus cause fiux-ratio anomali es 
(IChen et al.1 '2003: Wambsganss et al. 2005; Metcah'20Q5allbl: 
iMiranda k Macci o 2007; Puchwein k Hilbert 2009). 

In particular, IChen et al.l (|2003l l used the cross-section 
(optical depth) method to calculate the effect of both 
subhaloes intrinsic to the main lens and line-of-sight 
haloes. They found that the former would dominate the 
total lensing cross-section, although the exact percentage 
was highly sensitive to the spatial distribution of sub- 



structurefl; the latter - line-of-sight haloes, modelled as 
singular isothermal spheres - would c o ntribut e to ^ 10% 
of the total perturbation. iMetcalj (|2005al l performed 
ray-tracing simulations for the line-of-sight lens population 
(lO^M© ^ m ^ lO^Mo) in a ACDM universe, which he 
compared to several observed systems with measured cusp- 
caustic ratios. Assuming that haloes have Navarro , Frenk k 
White (NEW) profiles (Navarro et allll996l . Il997^. Metcalf 
found that the predicted abundance of line-of-sight haloes 
was enough to explain the observed flux-ratio anomalies 
of several representative cases. With a slightly different 
approac h and assuming hal o es to be singular isothermal 
spheres, iMiranda k MacciEl (|2007l ) found that, with the 
additional contribution from haloes along the line of sight, 
the observed flux-ratio anomalies can be reproduced with a 
high confidence level. 

In this work, we re-examine the perturbing effect of 
haloes along the entire line of sight from the source to the 
observer by using A^-body simulations to generate strong- 
lensing sight lines and quantify the flux-ratio distributions 
for multiply- imaged sources. In §2, we introduce our method 
for tracing light deflection through multiple lens planes. 
In §3, we present a summary of the cusp-caustic viola- 
tions arising from simple perturbation scenarios (varying 
the density profiles, angular positions and redshifts of the 
perturbers). In §4, we describe our method for generating 
"lensing lig hcones" in the Millennium II ACDM A^-body 
simulation (|Bov1 an-Kolchin et "aD l2009l l. Results from the 
analysis of these lensing cones are given in §5. In §6 we 
present results using a Monte- Carlo approach to account 
for haloes below the mass resolution limit of Millennium 
II (^ 10^ Mq). Our conclusions are given in §7. The 
cosmology of our lensing simulations is the same as that 
used for the Millennium-II simulation, with a matter density 
Qm = 0.25, cosmological constant — 0.75, Hubble con- 
stant h — ifo/(100 kms~^ Mpc~^) = 0.73 and linear fiuctu- 
ation amplitude as = 0.9. These values are consistent with 
cosmological constraints from the WMAP 1 - and 5-year data 
analyses (Spergel et al. 2003: [Komatsu et^a l. 2009), but dif- 
fer at about th e 2a level from more recent WMAP 7-year 
determinations (jKomatsu et al.|[201lh . This small offset is of 
no consequence for the topics addressed in this paper. 



2 SIMULATIONS OF LIGHT DEFLECTION 
THROUGH MULTIPLE LENS PLANES 



In Chapter 9 of ISchneider et al.1 (|l992l l. the authors present 
the theory of light deflection through multiple lens planes. 
As shown in Eig. [U haloes are projected near the line of 
sight at all redshifts between the observer and a source at 
redshift Zg. The angular position of the source is denoted 
by ^ and its final image position is denoted by 0* . The 
haloes are assumed to be distributed in N lens planes, each 
at redshift Zi (i=l,2...A/', and Zs — ^^a^+i). As the light ray 
passes through each plane, the image position OiJ^i (where 



^ Apart from quasar images' flux-ratio anomalies, another 
promising method is to use surface brightness anomalies of lensed 
galaxies to identify substructures and constrain their properties, 
see e.g., Vegetti k Koopmansl (120091) : IVegetti et al.l (l20ia V 



A similar conclusion was reached bv lXu et al.l (I2OO9' '2010^: see 
also Nieren berg et al.l (|201lh for the observed spatial distribution 
of luminous satellites in early- type galaxies. 
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primary lens Source 




Images of any given background source can be accu- 
rately and efficiently identified using the Newton-Raphson 
method, once the mapping relation (the overall Jacobian 
matrix As — is obtained. We determine the Jacobian 

matrix numerically, as follows. For each lens plane, a rigid 
grid of 1000 x 1000 is applied to cover a central region of 
5" X 5". Source positions I3n{0i) that correspond to grid 
points Oi in the first lens plane (which is also the final image 
plane) are calculated through the multi-plane lens equation 
(Eq. [1]). An arbitrary light ray propagating through a lens 
plane will not necessarily hit a grid point of the mesh that 
covers that plane. Therefore, for any given position Oi of 
the i-th lens plane, the defiection angle ai{Oi) is obtained 
by linear interpolation using the values for the four nearest 
grid points. Once I3n{0i) is obtained, the Jacobian matrix 
^s(^i)(= d/SN/dOi) for the image plane grid points Oi can 
be derived using finite differencing with the five-point stencil 
method (accurate to the third order). Again, linear interpo- 
lation is used to obtain each element of the Jacobian matrix, 
As{0), at any given image position 0. The corresponding im- 
age magnification /i(^) is then given by /i = detA~^. 

We use our multi-plane ray tracing code with a resolu- 
tion of 0.005^^ per pixel in the lens and image planes, which 
we find sufficient to accurately reproduce the lensing prop- 
erties of a number of simple analytical cases. 



Figure 1. An illustration (top) of a light ray propagating through 
intergalactic space from a lensed quasar to an observer. The pri- 
mary galaxy/halo at intermediate redshift causes image splitting 
due to its strong lensing effect. Both intrinsic substructures (satel- 
lite subhaloes and galaxies) in the primary lens and intergalactic 
haloes along the line of sight perturb the lensing potential and 
give rise to anomalous flux ratios between the images. The corre- 
sponding illustration of ray tracing through multiple lens planes 
is given in the bottom panel. 



ANOMALOUS FLUX RATIOS AND 
CUSP-CAUSTIC VIOLATIONS 



The cusp-caustic relation (iBlandforc 


k Naravan Il986l: 


ISchneider k Weisslll992l:IZakharovll 19951: 


Keeton et al.ll2003ll 



is defined as: 



Rcu 



(3) 



the light ray intercepts the plane) in the (i+l)-th lens plane, 
which is also the source position I3i for the i-th plane, is 
related to I3~^ and O'^ by the lens equation: 



= ft = 1 



3 = 1 



Di 



(1) 



is the defiection 



where On+i — I5n — I^^ ^ and Oi — 0^ . a 
angle a light ray undergoes in the j-ih plane at Oj. D^+i and 
Dj^i^i are angular diameter distances between the (z+ l)-th 
plane and the observer, and between the {i + l)-th plane 
and the j-ih plane, respectively. Dn+i = Dg is the angular 
diameter distance of the source. The Jacobian matrix Ai of 
the mapping between Oi and the source position /3i for the 
i-th plane is given by: 



A^. 



dOl dOl 



E 



■[ Di+i Q0. Q0^ 



(2) 



and An = As — is the overall Jacobian matrix, describ- 
ing the mapping relation between and 0* . 



when /itotai = |ma| + |mb| + |mc| oo. /i denote the mag- 
nifications of the three closest images (A, B and C) of a 
background point source located near a cusp of the tan- 
gential caustic (as shown in Fig. [2]). Observationally, source 
positions cannot be directly measured - instead, an image 
opening angle is often used as an indicator of the proxim- 
ity of a source to the nearest cusp of the tangential caustic. 
This opening angle is measured between lines joining the 
centre of the lens to the two outer images A and C. As the 
source moves outwards (towards the nearest cusp), 0, 
Mtotai oo, and i?cusp will go to zero asymptotically. This 
relationship holds for any smooth lens potential. 

Fig. [2] illustrates how perturbing structures change the 
cusp-caustic relation. The upper panels show the critical 
curves in the image plane, and the bottom panels are con- 
tour maps of i?cusp for sources within the tangential caustic 
in the source plane. Left and right columns show smooth 
lens potentials and lens potentials with substructures, re- 
spectively. Substructures located near the critical curve will 
affect images nearby and result in significantly larger values 
of i?cusp, violating the predicted ratios of image magnifica- 
tions (fiuxes) given by Eq. (|3]). 
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Table 1. Four-image quasar lensing systems with ^ 120° 
measured at radio wavelengths (CLASS). 



smooth 



T 

substructure 




Figure 2. An illustration of how the presence of substructures af- 
fects the cusp-caustic relation. The upper panels show the critical 
curves in the image plane; the bottom panels are the contour maps 
of i^cusp for sources within the tangential caustic in the source 
plane. Squares indicate positions of close triple images and the 
corresponding sources in the two planes. The image opening angle 
is labelled for one case in the top left panel. The left column 
shows cases with a smooth lens potential. In the right column, 
we show cases where substructures are present. The cusp-caustic 
relation is violated when a perturbing structure is projected near 
the image positions around the critical curve (see text). 



3.1 Observational samples 

Multiple images of lensed quasars with small are ideal 
cases to examine violations of the cusp-caustic relation and 
can be used to put constraints on the properties of perturb- 
ing structures. This is especially true when their fluxes are 
measured in the radio and mid-infrared, as the interpreta- 
tion of optical and near-infrared flux ratios is complicated 
by stellar microlensing and dust extinction. At the present 
time, only five cusp-geometry lensing systems with image 
opening angle A^ ^ 90° are known. These were used for 
statist ical co mpariso ns to the simulations in our previous 
work (|Xu et al. 2009. |201Q| ). All five cases have surprisingly 
large i?cusp values which are difficult to explain with simple 
smooth lens models. Of these five (flux-ratio measurements), 
two that were obtained in the optical have been proven to 
be affected by microlensing; the other three were from the 
CLASS survey (Browne et al. 2003; Mvers et al. 2003) at ra- 
dio wavelengths and are thought to be more secure cases of 
perturbations due to substructure s in the lens. 

Table 2 of IChen et al] (|201lh lists all of the currently 
observed Rcusp-^0 pairs for systems with four distinct point- 
like images of quasars lensed by one single galaxy. Our Ta- 
ble [1] lists those with their flux ratios measured in the radio 
and image opening angles A^ ^ 120°. We compare our sim- 
ulations with this observational sample of lenses, with no 
additional selection criteria. 



Systems 



AO (o) Rc 



Reference 



(1) B2045+265 

(2) B0712+472 

(3) B1422+231 

(4) MG0414+053 

(5) B1555+375 



34.9 0.501 



76.9 



0.255 



77.0 0.187 



101.5 0.227 



102.6 0.417 



Fassnacht et al. (1999) 
J<oogrnari^^l^ (^003 ) 
Jackson et al. (1998) 
J<oogrnari^^l^ (^003 ) 
Patnai k et al. (1999) 
Koopman^^l^ (^003) 
Hewitt et al. (199^ 
Katz et al. (1997) 
Marlow et al. (1999) 
J<oayman^^l^ (^003 ) 



Notice that four of the five lenses listed in Ta- 
ble 1, namely, B0712+472, B1422+231, B2045+265 and 
MG0414+053, are reported as having visible companions 
(satellite s /group galaxies) projected near the main lensing 



(satellite s /group galaxies) projected near tne mam iensm^ 
galaxies jP assnacht Lubinll200l lHogg Blandfordl[l99l 
Grant et a l. 2004; McK ean et all l2007l : iFalco et al.1 Il997l l. 



Fitting the observed image positions using a singular isother- 
mal ellipsoidal model yields velocity dispersions a ranging 
from 200 ~ 400 k m/s, and axis ratios q between 0.75 — 0.9 
(|Sluse et al.l [201lV exce pt for B1555, wh ich requires a = 133 
km/s and q = 0.45 (Marlow et al.lll999l ). 



3.2 Statistical measures for the cusp-caustic 

violation: P(^ i?cusp| A6' ± 2.5°) and P^°(i?°usp^) 

Given a simulated lensing system, we compare to the obser- 
vations in Table 1 by generating a large number of realisa- 
tions of background sources with ^ 120°. We calculate 
^cusp for each realisation and evaluate P(^ i?cusp| A^±2.5°) 
for this ensemble of realisations. This is defined as the prob- 
ability for i?cusp, measured for sources with image opening 
angles G [AO - 2.5°, AO + 2.5°] (i.e. within a five-degree 
opening-angle span centred at A^) to be larger than a par- 
ticular threshold value. Lenses with more perturbations will 
result in large i?cusp values for many source positions and 
thus have a higher P(^ i?cusp|A^ ± 2.5°) than lenses with 
fewer perturbations. 

We illustrate our use of the P(^ i?cusp|A^ ±2.5°) in 
Fig. [3l The top panel shows a typical example of a close 
triple image configuration for cusp sources with A^ ^ 120°. 
In this case, the lensing galaxy has a (smooth) singular 
isothermal ellipsoidal (SIE) profile (seel Keeton KochanelJ 
(1998) for notations of 6/, 6sie, ^3, and so herebelow) with 
lensing strength bi = 0.6^^ and axis ratio ^3 = 0.8, and is 
located at redshift Zd = 0.6; the source redshift is = 2. 
The corresponding contour map of P(^ i?cusp|A^ ± 2.5°) 
in the i?cusp-A^ plane is given in the middle panel. 
Also plotted are the radio measurements for the cur- 
rently best available sample (listed in Table [1]). These are 
clearly inconsistent with the smooth-lens i?cusp distribution. 

When we include the substructures within the lens- 
ing galaxy and its dark matter halo, the regular i?cusp 
distribution for a smooth lens potential disappears. The 
bottom panel in Fig. [3] shows the average distribution of 
P(^ i?cusp| A^zb2.5°) when including the subhalo population 



Line-of-sight lensing and flux anomalies in LCDM 5 




from the Aquarius simulations (IXu et alJ l2009l ) . At small 
A^, violations are more significant than on larger scales. The 
smallest i?cusp measured among all observ e d cus p-caustic 
systems is 0.187 (from B1422). In IXu et al.1 (|2009l ). we cal- 
culated P^^{R[ 



.0.187^ 

cusp 



Figure 3. Top panel: Close triple image configurations for a SIE 
lens with 6/ = 0.6^^ and qs = 0.8. The critical curve of the lens 
is shown in black and its caustics in grey. The regions sampled 
by "cusp sources" are shown in red, and the corresponding dis- 
tributions of the three images are shown as green, pink and blue 
regions around the critical curve. Middle panel: the correspond- 
ing probability contour map for Rcusp\^0 ± 2.5°). Contour 
levels of 1%, 2%, 5%, 10%, 20% and 50% (from top to bottom) 
are plotted. Blue squares are the five radio measurements (with 
^ 120°) so far available. Bottom panel: the average probabil- 
ity contour map (for A^ ^ 90°) of violatio ns due to s ubstructures 
in the Aquarius haloes (using results from lXu et al.|[2009. ). 



which is the probability for i^cusp to be 
larger than or equal to 0.187, computed over all realizations 
with AO ^ 90°. P^°(i?cusp^) was found to be -10%. We 
concluded that it is difficult to explain the observed i?cusp 
distribution (especially at larger A^) with a subhalo popu- 
lation similar to that produced in the Aquarius simulations. 
This motivates the search for other sources of perturbations 
to the lens potential. 

In this work, we use P^^{Rcus^) as an overall estimate 
for the probability of observing cusp-caustic violations, in 
order to compare with our previo us work. The value of 90° 
is chosen in order to compare with lAmara et al] (|2006l ) , who 
adopted A^ ^ 90° to select cusp-like lenses which can be 
best used to test the cusp-caustic violation. Varying the up- 
per limit of A^ will change the probability of cases with 
i^cusp ^ 0.187. However, our final conclusion (in §7) is based 
on a statistical argument that only applies to the observed 
lenses with AO ^ 90°. 



3.3 Simple perturbation scenarios with different 
halo redshifts, masses, profiles and 
concentrations 

A number of parameters determine the importance of these 
perturbers for creating violations to the cusp-caustic rela- 
tion: most significant are their masses, density profiles, red- 
shifts and impact parameters to the line of sight. Before 
presenting the results from general lines of sight taken from 
A/'-body simulations, we first show several simple perturba- 
tion scenarios to illustrate, individually, the eflFects of these 
difl[erent parameters, in the case of a single perturbing halo. 

Fig. 3] shows critical curves and caustics produced by a 
main- lens potential of an isothermal ellipsoid with bi = 0.6^^ 
^3 = 0.8 and core size so = 0.05^^ located at Zd = 0.6 for a 
source at redshift Zs = 2, plus a perturber of m = 10^° M© 
modelled with a truncated singular isothermal sphere. The 
panels in this figure correspond to diflFerent scenarios. In the 
upper row the perturber's angular position is fixed (out- 
side the tangential critical curve), and we vary its redshift: 
z = 0.4 (foreground), z = 0.6 (in the main- lens plane) and 
z — lA (background). In the lower row, we fix the redshift 
of the perturber to the main lens plane [z — 0.6) and change 
its impact parameter, such that it is projected within, on top 
of and outside the tangential critical curve (left to right pan- 
els, respectively). Wiggles and swallow tails are introduced 
to the critical lines and caustics by the added perturbing 
structure; massive perturbers (or those with compact den- 
sity profiles) can even cause a secondary set of criticals and 
caustics. Images located around these wiggle features violate 
the cusp-caustic relation most strongly. 

Fig.[5]shows the contour maps of P(^ i^cuspl A^ zb 2.5°) 
in the Rcms^-AO plane for the diflFerent scenarios above. Vi- 
olation patterns as a function of image opening angle A^ 
vary with the positions and redshifts of the perturbers. No- 
tice that at redshifts greater than that of the primary lens, 
the cone of light rays starts decreasing in size towards the 
source. This means a "background" perturber that appears 
to be projected close to the critical curve (where images nor- 
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mally form) could actually be far away from the light ray. 
Such perturbers would be less effective in causing conver- 
gence fluctuations than their foreground counterparts (they 
would still contribute to the shear field). However, depend- 
ing on the distribution of the primary lens and the source, 
there could be many more background structures affecting 
the light ray than those in the foreground (see §5). 

The mass and density profile of a perturber also af- 
fect the production of flux-ratio anomalies (and cusp-caustic 
violations) by altering the effective cross-section. Singu- 
lar isothermal spheres have been found to be a good 
approximation for the inner density profiles of relatively ' 



approximation tor tne inner density pronles oi relatively 
massive haloes (iRusin et al.ll2003l: iTreu KooDmaii3l2004l : 



iRusin k Kochanekl l2005l : iKoopmans et al.1 l20Q6l ). where 
baryons are thought to dominate their central potentials. 
This effect may be less important in smaller haloes, where 
the density profile is more likely to be well approximated 
by the NFW distribution c haracteristic o f CDM haloes 
in A^-body simulations (Navarro et al.|[l997l l. Nevertheless, 
there is still much controversy whether observed low-mass 
haloes around d warf galaxies have core- like shallow profiles 
(lOh et al.ll2QlTI ). 

In the simple scenarios presented below and in our 
line-of-sight lensing simulations (see §5 and §6), we model 
perturbing haloes either as truncated singular isothermal 
sphere^ or truncated NFW profiles, normalized with their 
masses and truncated at their virial radii. We follow the 
convention of defining the virial radius as r2oo, the radius 
within which the mean halo density is 200 times the critical 
density of the Universe (at the appropriate redshift z). The 
mass enclosed within r2oo is denoted as M200 • For the NFW 
profile, the concentration parameter is C200 = T2oo/ts, where 
Ts is the scale radius. This parameter is thought to correlate 
with mass M200 and redshift z. A number of concentration- 
mass relations have been proposed in the literature, based 
on A/'-body simulations. 

In this work, we a dopt the concentration-mass relation 
of iMaccio et al.l (|2008l l (hereafter M08) wherever we model 
perturbers as truncated NFW profiles. The fitting formula 
(for a WMAP-1 cosmology, close to that of the Millennium- 
II simulation) is given by: 



C2oo{M2oo,z) 



10" 



Mo 



(4) 



where H''{z) = Hg[QA + ^^m(l + z)^]. 

The concentration-mass relation of lBullock etHI (|200lh 
was used by iMetcalj (|2005al lbh to study how line-of-sight 
haloes (lO^M© ^ m ^ 10^ M©) contribute to the flux 
ano maly problem. The adopted fitting formula was given 
by (|Metcaljl2005bh : 



C200(M200, z) - 
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M2, 



(5) 



1 + z VlQi^M©^ 

To compare with iMetcalj (|2005al lbh. we also perform our 
analysis using this alternative concentration-mass relation, 
hereafter referred to as B01-M05. 

Fig. [6] shows how these different assumptions for the 

^ A singular isothermal sphere may not be a realistic model for 
small haloes. We adopt this model for ease of comparison with 
previous work, e.g. iChen et al.l (|2003|). 



mass, density profile and concentration-mass relation of 
a perturber change the total critical curves produced by 
a primary lens at z = 0.6 and the perturbing halo lo- 
cated at the same redshift and with a mass of M200 = 
[lO^Mo, lO^Mo, lO^^M©]. A different density profile (a 
truncated singular isothermal sphere, a truncated NFW pro- 
file with the M08 concentration-mass relation and a trun- 
cated NFW profile with the B01-M05 concentration-mass 
relation) is assumed in each row of Fig. [S] Different distor- 
tions to the critical curve correspond to different levels of 
violations in the smooth- lens flux-ratio relationship. 

The mass dependence of the violation pattern has 
been studied systematically, with results presented in §6, 
which also includes a discussion of effects from different 
concentration-mass relations and from allowing scatter in 
the concentration on the overall cusp- violation probabilities. 

In this section we have illustrated the effects of vary- 
ing the redshift, impact parameter, mass and density profile 
of a single perturbing halo. In practice, perturbations could 
arise anywhere along the line of sight and from many differ- 
ent haloes. The overall perturbation is far more complicated 
than any of the simple cases presented here. In the following 
sections, we use cosmological A/'-body simulations to obtain 
self- consistent and realistic distributions of perturbers along 
strong lensing sight lines, and estimate the net perturbation 
and the likelihood of the observed flux ratio violations. 



4 LENSING LINES OF SIGHT FROM 
COSMOLOGICAL SIMULATIONS 

4.1 Constructing lensing cones from MS-II 

The Millennium-II simulation f MS-II jBovlan-Kolchin et al.l 
2009) is an A^-body simulation of a cubic cosmological 
volume with a comoving side length of 100h~^ Mpc, at 
a spatial resolution of lh~^ kpc and mass resolution of 
6.89 X lO^/i~^M0. The cosmological parameters of MS-II 
are the same as those of the earlier Millennium and Aquar- 
ius simulations, consistent with the WMAP-1 results. MS-II 
provides us with the large-scale distributions of a cosmologi- 
cal sample of dark matter haloes. When tracing lensing sight 
lines through this simulation, we use the following method to 
determine where haloes cross the past light cone of a fiducial 
observer (for more details, see ^Angulo 2008). 

We start by replicating the 100/i~^ Mpc simulation box 
in its X, Y and Z dimensions, as many times as we need 
to cover the desired redshift range (along the sight line) and 
angular size (transverse to the sight line) . For computational 
efficiency we only let the combined box go to the source 
redshift in the X and Y dimensions, and keep the number of 
replications in the Z dimension to a minimum. As illustrated 
in Fig. [71 the observer is located at the origin (0, 0, 0) in one 
corner of this replicated box. Assuming a source redshift of 
= 2, the total dimension of the combined box is chosen to 
be 38 X 38 X 8, in units of one Millennium-II simulation box. 
The position angles {0, cj)) of a given line-of-sight vector are 
defined as the angles measured from the ZX- and Xy-plane, 
respectively. The simulated sky into which we trace sight 
lines is then two 2° x 30° stripes, which cover 10° ^ 6* ^ 40° 
and 50° ^ 6* ^ 80°, and 10° ^ ^ 12°. Directions along the 
X and Y axes (with ^ 10°) and along 40° ^ 6* ^ 50° have 
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Figure 4. Critical curves (six panels on the left) and caustics (six panels on the right) produced by the same smooth lens potential 
plus a m = lO^'^ Mq perturber (indicated by red squares), modelled by a truncated singular isothermal sphere. Top row: the perturber's 
angular position is fixed and its redshift set at z = 0.4 (in the foreground, left panel), z = 0.6 (in the primary lens plane, middle panel) 
and z = 1.4 (in the background, right panel). Bottom row: The redshift of the perturber is fixed (z = 0.6) and its impact parameter is 
changed from inside the tangential critical curve (left panel, labelled as "In"), to overlapping (middle panel, labelled as "On") to outside 
(right panel, labelled as "Out"). Parameters are noted in each panel. Note the wiggles induced in the critical curves and the production 
of secondary critical curves and caustics. 
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Figure 5. Corresponding probability contour maps of P(^ i?cusp|A^ ± 2.5°) for cases presented in Fig. 31 Symbols and contour levels 
are the same as in Fig. [S] The top row presents cases where a perturber (of 10-^^ Mq) is projected at the same angular position as shown 
in the upper panels of Fig. [H but located at different redshifts: z = OA (foreground), z = 0.6 (in the main-lens plane) and z = lA 
(background); the second row shows cases where the perturber is located at z = 0.6 but projected at three different angular positions as 
shown in the lower panels of Fig. [Jl 
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Figure 6. Critical curves for smooth lens potential with a single perturbing halo (located at z = 0.6), plotted as a red square. In 
columns from left to right the perturber has a mass of lO^M©, lO^M©, and W^^Mq, respectively. The rows from top to bottom 
correspond to different assumptions for the density profile: a truncated singular isothermal sphere, a truncated NFW profile with the 
M08 concentration-mass relation, and a truncated NFW profile with the B01-M05 concentration-mass relation, respectively. 



been excluded to avoid significant repetition of structures in 
the replicated box. 

Haloes at each simulation snapshot (corresponding to 
a particular redshi ft) are identified using the Friends-of- 
Friends algorithm (|Davis et al.l Il985l l . Haloes also contain 
many subha loes; these are identified using the subfind 
algorithm ( Sprin gelll2005[ ) . The minimum mass of subhaloes 
resolved by the simulation is 1.4 x lO^/i^^M© (correspond- 
ing to 20 particles) . Haloes at different snapshots are linked 
together by an algorithm for defining their merging history 
(|Hellv et al.L 2QQ3). We follow haloes in each of these merger 
trees and predict their trajectories (in the replicated box) 
between every two adjacent snapshots. In this way, we can 
find the exact redshift and comoving position of a halo at 
the moment it crosses the past light cone of the observer. 
When a halo crosses the light cone, all its subhaloes are 
assumed to cross at the same redshift. We assume that 
the relative positions of these subhaloes at the light-cone 
crossing time are the same as in the previous snapshot. 



Hereafter, we will use the term "lensing cone" to refer to 
the observer's light cone that encloses a particular lensing 
sight line towards a certain direction in the sky (and out 
to the source redshift). All haloes that cross the past light 
cone are checked to see if they are physically within a given 
lensing cone. If so, their positions, redshifts, masses and half 
mass radii are stored for lens modelling. All lensing cones are 
50^^ X 50^^-wide, out to a source redshift Zs = 2, and contain 
a primary lens around redshift Zd = 0.6 (typical source and 
lens redshifts for the observed quasar lensing systems). 

To build up our lensing cone catalogue, we randomly 
select about 300 directions in the 2° x60° simulated sky, each 
of which goes through at least one galaxy-scale halo with a 
mass above 10^^ Mq located at redshift |z — 0.6| ^ 0.02 in 
the replicated box. This ensures that the primary haloes we 
select are responsible for producing multiple images of the 
= 2 background sources. We have confirmed that these 
^ 300 randomly selected primary lenses are representative 
of mass and circular velocity distributions of ^23,000 haloes 
that meet the same selection criteria in the simulated sky. 
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Figure 7. The geometry of the rephcated box for hght-cone gen- 
eration: the MS-II simulation box of 100h~^ Mpc is repeated in 
its X, Y and Z dimensions as many times as needed to cover 
the desired redshift range and angular size. For a source redshift 
of Zs = 2, the total dimension of the combined box is set to be 
38 X 38 X 8, in units of one Millennium-II simulation box. An ob- 
server is put at the origin (0, 0, 0) of this box. The position angle 
pair (6, (f)) of a given line-of-sight vector are defined as the angles 
measured from the ZX- and Xy-plane, respectively. The simu- 
lated sky we have looked at is then two 2° x 30° stripes, which 
cover 10° ^ 6* ^ 40° and 50° ^ 6* ^ 80°, and 10° ^ (/> ^ 12°. 
Directions along the axes (with ^ 10°) and along ~ 45° have 
been excluded to avoid significant structure repetition. 



As stated in §4.2, velocity dispersions of the simulated 
main lensing haloes range from 200kms~^ to 300kms~\ 
comparable to our sample of observed lenses (see the end 
of §3.1). As the selection function of the observed lenses is 
hard to define, we assume that they are a random sample 
of haloes in the ranges of velocity dispersion and redshift 
given above. Furthermore, by imposing a lower mass limit 
of lO^^/i-^Mo on the main lens, we have excluded cases 
where two or more less massive haloes aligned along the 
same line of sight produce a comparable strong-lensing sig- 
nal. However, the chance is small for two foreground haloes, 
both more massive than lO^°/i~^M0, to be well aligned to 
jointly lens a background galaxy. Emprically, such "three- 
dimensional" lenses are rare: in the CLASS survey, only one 
probable su ch case has been reported out of a t otal of 22 
candidates (Cha e et al.ll200ll : lAugusto et al.ll200lh . 

Fig. [8] shows the geometry and halo distribution of an 
example lensing cone. All haloes within a given cone are 
used for the lensing calculation. On average each lensing 
cone (of 50'' x 50'') contains about 10,000 (12,000) haloes 
(subhaloes). Within a projected central region of i? ^ 5" 
for strong lensing, there are on average ^ 300 haloes with 
m > 10^ h~^MQ directly contributing to the convergence 
field. The rest are distributed further out (in projection) 
and contribute to the shear field of this region in the same 
way as point masses. 

4.2 Ray-tracing through MS-II lensing cones 
and line-of-sight lens modelling 

To carry out calculations for multi-plane light deflection, 
we assume 60 lens planes distributed with equal spacing in 
redshift between the observer and the source at = 2. In 
each of these lens planes, a region of 5" x 5" around the 
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Figure 8. The halo distribution within an example lensing cone 
m a slice of depth 3800/1" ^ Mpc (in co-moving distance), out to a 
redshift of 2; = 2. The cone covers a region of 50^^ x 50^^ In all four 
panels, haloes and subhaloes are plotted as black squares; those 
more massive than 10^^ Mq are shown by red squares. The central 
region of radius R ^ 5^^ is indicated with blue lines in both top left 
and bottom right panels. Top left: the projected cone geometry 
in comoving distance. Top right: same lensing cone, shown in 
the redshift — physical distance plane. Bottom left: the central 
4// ^ 4// j-ggiQn of the light cone projected in the sky, where the 
main tangential critical curves form; subhaloes are indicated by 
purple squares in this panel, circles in solid red and dashed black 
indicate the half-mass radii of haloes (and subhaloes). Bottom 
right: an expanded view, showing haloes and subhaloes projected 
within the central 40^^ x 40^^ of the light cone. 



line centre is covered by a 1000 x 1000 rigid grid in order to 
calculate the Jacobian matrix As (Eq. [2]) between the source 
plane and the final image plane. 

Haloes within a lensing cone are projected into these 
lens planes according to their redshift s. The main lens halo 
is modelled as an isothermal ellipsoid, for which a universal 
axis ratio (^3 = 0.8) and core radius (so = 0.05") are as- 
sumed. The orientation of the ellipsoid is randomly chosen 
in the interval of [0, 27r]. The lensing strength bsm (related 
with bi through h i — bsme/s in~^e, where e = (1 — ^3)"*^^^, 
see iKeeton Kochaneklll998h is derived through an empir- 
ical relationship between halo's virial velocity V200 and the 
veloc ity dispersion cts ie of the equivalent isothermal ellip- 
soid tehae et ahlboOGl l: 

"1 0.22 

(171 km/s ^ 1^200 ^ 563 km/s) 



<7SIE 



200 km/s 



1.17 V2i 



200 km/s 



(6) 

and bsm = 47r(crsiE/c)^i^ds/^s, where c is the speed of 
light and Dds and Ds are the angular diameter distances 
between the main lens and the source, and the source and 
the observer, respectively. The virial velocity V200 is ob- 
tained from halo mass M200 and its virial radius r2oo through 
V200 = G'M2oo/t2oo- Our requirement that the main lens 



be more massive than 10^^ h' 
from 200 km s"^ to 300 km s" 



''^Mq results in ctsie ranging 
^, with a weighted mean (by 
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cross-section, oc a^) of 222kms~^ corresponding to a lens- 
ing strength of b^m = 0.84^^ for our adopted lens and source 
redshifts {zd = 0.6, Zs = 2.0). 

Within each lensing cone, haloes with projected pro- 
files that are completely outside the central ^" x ^" region 
are treated as point masses. Those within this region are 
assigned a density profile: as described above, we investi- 
gate three distinct choices of this profile (singular isother- 
mal sphere, NFW with the M08 concentration-mass relation, 
and NFW with the B01-M05 concentration-mass relation). 
All halo profiles are normalized to their masses M200 and 
truncated at the virial radii r2oo ; subhaloes are truncated at 
two times their half mass radii. 

For each line of sight, deflection angles are individu- 
ally calculated for the equivalent isothermal ellipsoid of the 
main lens and for all line-of-sight (sub) haloes, and are tab- 
ulated to the meshes at different lens planes. Through ray 
tracing, source positions Pn that correspond to the final 
image plane Oi are identified, and the final Jacobian matrix 
As — dj^N /dOi is then derived using the finite differencing 
method. Images of a given source position are effectively 
found using the Newton- Raphson iteration method. 



5 RESULTS FROM THE MILLENNIUM II 
SIMULATION 

In order to calculate P(^ i?cusp|A^ ± 2.5°) - the probabil- 
ity distribution in the i?cusp — plane for individual lens- 
ing cones, we generate 10, 000 ^ 20, 000 cusp sources whose 
close triple images have image opening angles ^ 120°. 
We have also calculated P^°(i^cuip^) for cases with A6' ^ 90° 
as an overall estimate for cusp-caustic violations to compare 
with our previous work, in which only cases with A^ ^ 90° 
were examined for violations (caused by intrinsic substruc- 
tures within the main lens). 

To derive average violation probabilities over all sight 
lines, we weight P(^ i?cusp|A6' ± 2.5°) and P^°(i?cusf ) of 
each lensing cone by the quadruple-image lensing cross- 
section in the source plane (simply the fractional area within 
the tangential caustic) . We do not account for magnification 
bias among the cusp sources. 

Fig. [9] (upper panels) and Table [2] show that if perturb- 
ing structures are haloes (and subhaloes) distributed outside 
the main lens along the line of sight as in the Millennium-II 
simulation, where such haloes are resolved to lO^/i^^M©, 
they cause a non-negligible amount of cusp violations, com- 
parable to those due to the substructures in the lens itself. 
However, the violation pattern (as a function of A^) depends 
strongly on the density profiles applied to haloes projected 
near the centre of the line of sight. 



5.1 Effects from massive line-of-sight haloes 

We have also investigated the effect from line-of-sight haloes 
more massive than lO^^/i^^M©, which are most likely to re- 
tain a significant fraction of baryons in their dark matter 
potential wells. The chance of finding at least one of these 
massive secondary lenses intercepting a strong-lensing sight 
line (i.e. projected within a typical Einstein radius of l" 
around the line centre, out to a redshift of 2) is about 10%. 



Depending on their density profiles, compact haloes (e.g. sin- 
gular isothermal spheres) could generate severe astrometry 
anomalies with a probability of a few percent, while haloes 
with shallower inner profiles could not. 

Fig. [10] presents the peculiar image configurations for 
an example sight line. Four, six and eight images (excluding 
the central image) are produced when the source is located 
at different positions with respect to two sets of tangen- 
tial caustics. In this particular case, the second caustic is 
produced by a perturbing halo of 2 x 10^° M©, modelled 
as a truncated singular isothermal sphere, projected near 
the centre of the main lens. Such peculiar image astrome- 
try has already been proposed and used to constrai n density 
profi l es of massive inter gala ctic objects fe.g.. Wvithe et al.l 
I2OOII : IWilkinson et ahllioOlh . 

5.2 Substructures inside line-of-sight haloes 

To investigate the effect of substructures inside haloes along 
the line of sight, we have excluded all substructures from our 
Millennium-II lensing cones and calculated violations due to 
"smooth" line-of-sight haloes alone. Table [2] lists violation 
probabilities in this case (for different halo density profiles) . 
The relevance of subhaloes to lensing fiux-ratio anomalies 
strongly depends on their assumed density profiles. Com- 
parison of the first and second rows of Table 2 shows that 
for m > 10^ Mq, NFW-like substructures within hne-of- 
sight haloes are responsible for causing a few percent of the 
cusp-caustic violations. 



5.3 Background vs. foreground 

We have separated line-of-sight haloes that are distributed 
in front of and behind the main- lens plane (zd — 0.6). The 
violation probabilities of these two groups are listed in Table 
[21 calculated excluding their subhaloes. A higher violation 
probability is found caused by haloes in the background than 
in the foreground, as more haloes intercept the light rays be- 
hind the main lens plane, given a typical lensing geometry 
{zd = 0.6 and Zs = 2). It is interesting to notice that viola- 
tions from the foreground and the background roughly add 
up to the total violations due to haloes along the entire line 
of sight (second row of Table [2]). The ratio between viola- 
tions from the foreground and from the entire sight line is 
close to 2:5, which is the ratio between the comoving radial 
distances for Zd = 0.6 and = 2. 



6 RESULTS FROM MONTE-CARLO HALOES 
WITH A SHETH-TORMEN 
MASS FUNCTION 

The Millennium-II simulation has a limited mass resolution. 
To investigate the mass dependence of the violation pat- 
tern below the halo mass of lO^/i^^M©, we have used a 
Monte- Carlo method to generate intergalactic halo popula- 
tions with masses IO^H'^Mq ^ m < lO^^/i^^M© (see §7 for 
discussion on adopting 10^ Mq as the lower mass limit). 
Thes e haloes are drawn from the Sheth-Tormen mass func- 
tion (|Sheth k Tormenll2002l ) generated with the code pro- 
vided bv lReed et all (|200/t ). 

We have randomly generated 200 lensing cones out to 
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Table 2. ^^*^(^cusp^) from cosmological simulations. 



line-of-sight perturbers 


truncated singular isothermal sphere 


truncated NFW (M08) 


truncated NFW (B01-M05) 


haloes + subhaloes 


7.8% 


7.3% 


12.8% 


haloes only 


7.0% 


5.2% 


9.7% 


background haloes 


5.0% 


3.2% 


5.8% 


foreground haloes 


2.0% 


1.9% 


4.0% 



Note:_P^i?2(jgp^) ~ 10% was derived using only substructure populations (m > 10^ h ^Mq) from the Aquarius simulations (|Xu et al.l 
I2009l . l20ld ). Cases here are using line-of-sight structures (haloes and subhaloes) from the Millennium II simulation (m > 10^ h ^Mq); 
subhaloes from the main lensing halo have been excluded. 
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Figure 9. Contour maps of the violation probability Rcusp\AO ± 2.5°). Symbols and contour levels are the same as in Fig. O The 
upper panels are for MS-II lensing cones with all line-of-sight haloes and subhaloes (m > 10^h~^ Mq), and the bottom panels are for 
Monte-Carlo lensing cones with line-of-sight haloes that follow the Sheth-Tormen mass function (m ^ 10^h~^ Mq). Three columns from 
left to right correspond to different assumptions for the density profile: truncated singular isothermal spheres, truncated NFW profiles 
with M08 and B01-M05 concentration-mass relations, respectively. 



Zs — 2, each of which contains a main lensing halo modelled 
as an isothermal ellipsoid at redshift Zd = 0.6. The lensing 
strength bsm is fixed to be 0.84^^, the same as the mean 6sie 
of the main lenses in the selected sample of the Millennium- 
II lensing cones. The axis ratio — 0.8, core radius 5*0 = 
0.05^^ and an orientation angle of 0.257r are also taken to be 
the same for all main lenses. 

In each realization of the lensing cone, line-of-sight halo 
positions are randomly generated, with number densities as 
given by the Sheth-Tormen mass function at the redshifts of 
the 60 lens planes used for the Millennium-II lensing cones. 
Table [3] lists the mean surface number densities of projected 
haloes in different mass decades, averaged over 200 lens- 



Table 3. Mean surface number densities of projected haloes out 
to z = 2.0 per decade of mass, averaged over 200 Monte-Carlo 
lensing cones. 



[10^107] [10^10^] [10^,109] [109,10^0] ^ 10^0 (/z-iMo) 
414 50 6 0.7 0.1 (/arcsec^) 



ing conefl Haloes projected within the 50^^ x 50^^-cone are 

^ These numbers roughly follow a power-law mass function of 
dn(m) = m~^-^ dm, m being the mass of haloes. 







Figure 10. Critical curves and caustics of an example MS-II lensing cone. The remarkable wiggle feature on the tangential critical curve 
is caused by a 2 x 10^^ Mq halo along the line of sight, modelled as a truncated singular isothermal sphere. Green dots in red regions in 
the subpanels are example source positions. Red and blue dots in the main panels are the corresponding image positions; red for negative 
parities, blue for positive parities. Panels 1-4 present different image configurations for a source located within the caustic region where 
five images would be produced; panel 5-8 are for a caustic source with seven images; and panel 9 is for nine images. 



saved, and those projected within the central 5" x 5^^-region 
are modelled with truncated singular isothermal spheres 
and truncated NFW profiles (using both M08 and B01-M05 
concentration-mass relations). Those further out are mod- 
elled with point masses. Cusp-caustic violations were identi- 
fied in the same way as for the Millennium-II lensing cones. 



6.1 Mass dependence of the cusp-caustic violation 

The lower panels of Fig.[9]show i?cusp|A6' ± 2.5°) con- 
tour maps using Monte- Carlo realizations of the line-of-sight 
halo population, with a lower mass limit of 10^h~^ Mq. As in 
the MS-II lensing case, the frequency of cusp-caustic viola- 
tions depends strongly on the assumed halo density profiles. 
Tabled also presents the values of P^°(i?cusp^) when this 
lower mass limit is increased, so that the mass dependence 
of the cusp-caustic violations can be seen. 

Applying truncated singular isothermal spheres yields 
relatively larger contributions to cusp violations from more 



massive haloefl (m ^ 10^^^^h~^ Mq). When truncated 
NFW profiles are assumed, lower mass haloes would also 
cause a significant amount of violations. As can be seen from 
Table [4l when the lower mass limit of line-of-sight haloes 
decreases from lO^h'^MQ to lO^h'^MQ (lO'^ /i-^M©), the 
corresponding violation probabilities, P^'^(i?cusp^)? increase 
by -2% (1%), -6% (4%) and 12% (7%) when assuming 
truncated singular isothermal spheres and truncated NFW 
profiles with the M08 and the B01-M05 concentration-mass 
relations, respectively. 

Comparing Table [4] with Table [2l it can be seen that 
our Monte- Carlo results are similar to those obtained using 
the Millennium-II lensing cones in the mass range above 
10^h~^ Mq. This suggests that the clustering of haloes is not 

^ Xu et al ] (l2009l ) estimated the total lensing cross-section (Jcs oc 
^SIE ^ ^perturbers(^) fo^ singular isothermal lenses. In this case 
(Jcs oc m", and o; is a positive value, hence the total lensing cross- 
section will be biased towards massive haloes. 
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Table 4. Violation probabilities for Monte-Carlo lensing cones with main lens parameter 6siE = 0.84^^ 



truncated singular isothermal sphere 8.2% 7.7% 6.5% 4.8% 2.8% 

truncated NFW profile (M08) 12% 9.1% 5.6% 2.5% < 1% 

truncated NFW profile (B01-M05) 23% 18% 11% 5.5% 1.9% 



a dominant eflPect in the production of flux-ratio anomalies 
for galactic-scale lenses. 

6.2 Dependence on the Einstein radius 

When the lower mass cutoff for main lensing haloes in the 
Millennium-II lensing cones is reduced from 10^^ Mq to 
10^^ Mq, the mean lensing strength bsiE of the equiva- 
lent isothermal ellipsoids decreases from 0.84^^ to 0.62^^ (cor- 
responding to crsiE = 190kms~^ for Zd = 0.6 and = 2 
using EqEl). Table [5] presents the P^°(i?cusp^) values that 
result from four different 6sie for the main lenses in our 
Monte-Carlo lensing cones. In addition to bsm = 0.62^^ and 
0.84^^, we have calculated violations for an arbitrary bsm of 
1" (1.5^^), which is about the mean (largest) Einstein radius 
of the observed systems listed in Table [1] 

As can be seen clearly from Table O systems with larger 
Einstein radii have higher violation probabilities. This is ex- 
pected, because close triple images normally form around 
the tangential critical curve at about the Einstein radius. 
Comparing to the case of a small Einstein radius, a larger 
value of this radius results in a higher chance for the image 
triple (of a given opening angle A^) to be intercepted by 
line-of-sight perturbers. 

6.3 Halo concentrations and mass function 

As we have shown above, the cusp- violation probability de- 
pends strongly on our assumptions about halo concentra- 
tion. The concentration-mass relation derived by Bullock 
et al. (2001) has a simple functional form (including red- 
shift evolut i on) a nd has been widely used in the literature. 
IColm et al.1 (|2004l ) investigated concentration parameters for 
haloes of lO^/i^^M© ^ m ^ 10^ H'^Mq and found this rela- 
tionship to be a good fit. However, these early simulations 
of dark matter halo formation had relatively low numerical 

resolution and this can introduce systematic errors. 

Mo re recently, a numb er of authors in c ludin g lNeto et al.l 

(I2OO7I I. iGao et alJ (I2OO8I ). IZhao et al.1 (l2009l l and M08 
(whose results are used above), derived concentration-mass 
relations from high-resolution cosmological A/'-body simu- 
lations. These studies are restricted to haloes with m ^ 
lO^'^ Mq but they exhibit systematic differences from 
the concentrations obtained by B01-M05. For this reason 
we show lensing results using the M08 relation, extrapolat- 
ing to lower masses when required, but exploring how the 
results change when this relation is varied by factors of a 
few. 

The B01-M05 concentration-mass relation overesti- 
mates the concentration of small mass haloes inferred from 
the extrapolated M08 relation by factors of 3 ^ 4 at z = 0. 
Therefore, we expect the violation probability to be larger 
for the B01-M05 relation than for the M08 relation. The 
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Figure 11. The violation probability ^^°(i^2usp^) changes 
with the application of different concentration-mass relations. 
P^*-'(i?^;^gp^) is given by the Y-axis. The X-axis indicates the as- 
sumed halo concentration C (at any given mass) normalized to 
C(M08) - the concentration predicted by the M08 concentration- 
mass relation. The Sheth-Tormen mass function is used to gen- 
erate line-of-sight halo populations. Values of P'^^{B^^^^) at 
C/C(M08)=0.5, 1.0, 2.0 and 3.0 are plotted as the four black 
squares, which are connected by the black solid line. Assuming 
that violations grow linearly with the number of perturbers, the 
dash and dotted fines plotting (100 ± 20)% x and (100 ± 10)% X 
the P^'^(^cusp^) values as shown in the black solid line, resemble 
violation probabilities under (100±20)%x and (100±10)%x the 
Sheth-Tormen mass function, respectively. Red, green and blue 
squares present results (under the Sheth-Tormen mass function) 
when allowing for scatter of concentrations (for haloes of a given 
mass) in form of Gaussian distributions with the standard devia- 
tion being 0.1, 0.2 and 0.3 dex around mean concentration values 
predicted by the MOB concentration-mass relation. P^*^(i^cusp^) 
derived from using the B01-M05 concentration-mass relation is 
also given, indicated by the purple horizontal line. 



scatter in halo concentration also affects the final cusp- 
violation probability. 

Fig. [11] presents the values of P^°(i?cusp^) induced by 
line-of-sight haloes assuming the Sheth-Tormen mass func- 
tion. To allow for possible uncertainties in halo concentra- 
tion, we also show results for the case when the concen- 
trations inferred from the M08 concentration-mass relation 
are multiplied by factors of 0.5, 1.0. 2.0 and 3.0. Varying 
amounts of scatter in concentration (for haloes of a given 
mass) are modelled assuming a Gaussian distribution with 
mean value equal to the M08 concentrations. As may be 
seen, the violation probabilities depend strongly on halo 
concentrations. Higher concentrations result in higher cusp- 
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Table 5. ^^^(^cusp^) ^^^^ Monte-Carlo lensing cones with main lenses of different Einstein radii 6sie: violations due to line-of-sight 
perturbers more massive than l{)^h~^ Mq. 



truncated singular isothermal sphere 6.4% 8.2% 11% 15% 

truncated NFW profile (M08) 9.5% 12% 15% 20% 

truncated NFW profile (B01-M05) 20% 23% 28% 32% 



caustic violation probabilities. A larger scatter in concentra- 
tion will also increase the violation probability. 

The halo mass function (the number density of haloes 
per unit volume per decade i n ma ss) influences the cusp 
violation probability. iMetcalj (|2005bl ) found that flux-ratio 
anomalies caused by line-of-sight perturbers not only depend 
strongly on the radial density profile of the haloes (their con- 
centration), but also on the primordial matter power spec- 
trum on small scales. Miranda & Maccio (2007) suggested 
that flux-ratio anomalies could be used statistically as a 
test of the behaviour of the matter power spectrum on small 
scales. We do not explore these effects here but in Fig. [TT] 
we show the result of using (100 ± 10/20)% x the Sheth- 
Tormen mass function, and assuming that violations grow 
linearly with the number of perturbers. 



7 DISCUSSION AND CONCLUSIONS 

We have examined the effects of intergalactic cold dark mat- 
ter haloes on flux-ratio anomalies in multiply-lensed quasar 
images by ray-tracing along strong- lensing sight lines that 
are either taken from the Millennium II simulation (for 
haloes and subhaloes with m > 10^h~^ Mq), or generated 
using the Monte-Carlo method assuming a Sheth-Tormen 
mass function (for haloes with m ^ 10^h~^ Mq). 

We use P^^{Rcus^), the probability for the cusp-caustic 
relation, i?cusp, to be larger than or equal to 0.187 - the 
smallest value of i?cusp measured for cusp-caustic systems 
to date (for the quasar B1422) - over all realizations with 

^ 90°, as a statistical measure of the cusp-caustic vi- 
olation probability. We have found that the mean viola- 
tion probability from intervening haloes depends strongly 
on t heir density profiles. 

I Chen et al.l |2003) assumed singular isothermal spheres 
for line-of-sight haloes and find that they only contribute to 
^ 10% of the total perturbation. Assuming the same halo 
density profile, we find that the cusp-caustic violation prob- 
ability caused by line-of-sight haloes with m ^ lO^/i^^M© is 
comparable to that caused by intrinsic substructu res within 
the m ain lensing halo (P^°(i?cusf ) ~ 8 % vs. 10%,IXu et al. | 
1 20 id) , which is in good agreement with Miranda & Maccio 
(|2007l ). The different results between Chen et al. (2003) and 
ours can be attributed to the drawbacks of their cross- 
section method, which underestimates effects from mo re so- 
phisticated perturbation scenarios fsee lMetcaill2005al ). 

When we assume truncated NFW profiles for the line- 
of-sight haloes (m ^ lO^/i^^M©), the violation probabil- 
ity, P^°(i?cusp'^), increases to 23% if we adopt the B01-M05 
concentration-mass relation and to 12% if we adopt our pre- 
ferred relation, that by M08. These values are larger than 
that due to the intrinsic subhalo populations alone. 



A typical NFW profile has an Einstein radius 3 4 or- 
ders of magnitude smaller than a singular isothermal sphere 
with a same mass. However, NFW perturbers in the mass 
range from 10^h~^ Mq to lO^^^'^ Mq cause more cusp 
violations than their singular isothermal counterparts. This 
may be due to the fact that in this mass range, perturbation 
in magnification (ratios) is mainly from fluctuations in the 
local density field that do not change the image positions. 
When comparing an NFW with a singular isothermal sphere 
of the same mass, we notice that the surface density distri- 
bution of the former exceeds that of the latter from a radius 
of ^ 0.001r2oo outwards, which means the NFW profile is 
more effective in introducing fluctuation to the convergence 
and thus causing flux-ratio anomalies. 

On the other hand, the deflection angle of a perturber 
of m - lO^~^/i"^M0 is always small (< O.OOl'' for a singular 
isothermal sphere locating at z = 0.6), until the perturber is 
massive (m > 10^^ h~'^ Mq) and compact enough (a singular 
isothermal sphere) to have a deflection angle (> 0.01^^) that 
can shift a nearby image to a new position with a differejit 
magnification from the primary lens fsee iMetcalJ l2005bl ). 
This can explain the larger violation probabilities (as shown 
in Table 4), caused by singular isothermal perturbers of 
m ^ lO^"^°/i~^M0 than by their NFW counterparts which 
are less effective in causing flux- anomalies due to shifting 
image positions. 

Another issue co ncerns the finite- source effect. 
iMetcalf k Amaral (|2012l ) pointed out that biased results 
about substructures could be drawn due to the point source 
approximation, which is used in this work. 

The radio-emission regions of observ ed quasars are 
estim ated to be ^10 p arsecs in extent (|Andreani et al.l 
1999 : IWvithe et al.|[2002 ). corresponding to - 0.001'' for a 
source at Zs — 2.0. When the perturbing mass drops down 
below lO^/i~^M0, the corresponding effective perturbing 
area decreases to < 0.001'' in radius for the perturber at 
Zd — 0.6, becoming smaller than an image with /i ^ 10 — 20 
(around the tangential curve) of the radio emission region 
of a background quasar. As a result, the induced mag- 
nification fluctuation would be smeared out (within the 
image area), and thus no significant image flux anomaly 
would be observed at radio wavelengths (but could still be 
seen in the optical /near- infr ared, which comes from mu ch 
smaller physical regions. See iMoustakas k, Metcalj|2003l for 
spectroscopic gravitational lensing). This is why we do not 
consider the violation probability produced by perturbing 
haloes below 10^h~^ Mq. As can be seen from Table 4, 
even if we neglect contributions from perturbers below 
lO^/i" ^Mq, we still find - 10% cusp-violation probability 
from line-of-sight NFW-like perturbers adopting the MOB 
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concentration-mass relation. 

Several other points are worth noting. Firstly, the vio- 
lation probability depends, of course, on the concentration 
of the halo, and both large halo concentrations and a large 
scatter in concentration will result in higher violation prob- 
abilities. Thus, it may be possible to use the statistics of 
flux-ratio distributions (measured in the radio) from large 
samples of lensed quasars to constrain the density profiles 
of low-mass dark matter h a loes. 

Secondly, in Xu et al. ' ('2009') we found that the viola- 
tion probability is higher for systems with larger Einstein 
radii because the mass fraction in dark substructures in- 
creases with radius. Here, the probability of a system vio- 
lating the cusp-caustic relation is also seen to increase with 
the Einstein radius but for a different reason: close triple 
images (with a given opening angle A^) that form at larger 
radii are more likely to be intercepted by line-of-sight per- 
turbers. Adopting bsiE = 1.0^^ the mean Einstein radius for 
the observed sample, we find that the violation probability 
P^°(i?cusp^) increases to 15% if the M08 concentration-mass 
relation is adopted. We also note that if we use the BOl- 
M05 concentration- mass relat ion, the violatio n probabilities 
for the representative cases in iMetcalj (|2005al ) can be repro- 
duced. 

Third, the Rcusp-^0 distribution varies with the el- 
lipticity of the main lens in such a way that lenses with 
higher ellipticity (smaller axis ratio (73) have larger Rcusv , 
and thus potentially lar ger P^°(i?cusp^) (|Keeton et al.ll2003l : 
iMetcalf Amarall2012l ). In this work, we adopt qs = 0.8 for 
the simulated main lensing haloes. For four of the five lenses 
shown in Table 1 (including all three with image opening an- 
gle ^ 90°), the axis ratios (0.75-0.9 as inferred from the 
singular isothermal ellipsoid models, see §3.1) are close to 
the value we adopt (0.8) in our simulations. In the absence 
of perturbing structures, changing the axis ratio to ^3 = 0.7 
causes a negligible increase in P^^{Rcus^)- 

Fourth, the probability that a massive halo (m ^ 
10^^ Mq) intercepts a galaxy-scale strong- lensing sight 
line with an impact parameter of ^ l" from the main lens 
centre is about 10%. T hese halos can pertu rb image fluxes, 
surface brightnes s fe.g., | Vegetti et al.ll2010l l and even image 
astrometry fe.g., IWvithe et alT200lh . Halos with compact 
density profiles (e.g., singular isothermal spheres) could gen- 
erate extra image pairs locally with an image separation of 
0.01^^ ^ 0.1^^ resulting in more than four bright images of 
a background quasar. With upcoming lensing surveys, ob- 
servations of peculiar image configurations could put better 
constraints on th e density profiles of these massive haloes 
(|Orban de Xivrv k Marshall 2009). 

Finally, for masses above lO^/i^^M©, we find that 
halo clustering has only a minor effect on the flux-ratio 
anomalies for galaxy-scale lensing systems. For a typical 
lensing geometry (with Zd = 0.6, Zs = 2), the overall 
perturbation produced by background haloes (behind the 
main lens) is larger than that caused by foreground haloes. 

To summarise, in this work we have calculated the cusp- 
caustic violation probability, as measured by ^^^(i^cusp^), 
produced by line-of-sight dark matter haloes. The value 
of P^°(i?cusp^) strongly depends on halo density profiles, 
specifically on concentration, in the case of NEW per- 



turbers. When the c oncentration-mass relation proposed by 
iMaccio et aP (|2008l ) is used, the value of P^°(i?cusp'^) pro- 
duced by all line-of-sight perturbers is found to be ^ 10 — 
15% (correspondi ng to a mean Einstein radius of 0.8 — 1.0^^) 



In previous work (IXu et al.| 



simulations fepringel e t al 



2009 



I2OI0I ) using the Aquarius 
20081 ), we found that the con- 



tribution to P^°(i?cusp^) from substructures inside the main 
lensing haloes also amounts to P^°(i?cusp^) 10 — 15% (cor- 
responding to a mean Einstein radius in the same range as 
above). Summing up both contributions, the total violation 
probability could reach 20 — 30%. 

There are five observed cusp-geometry lensing systems 
whose triple images have opening angles A^ ^ 90° . Of these 
five, the three radio lensing cases (B0712+472, B1422+231 
and B2045+265 as listed in Table 1) show firm evidence 
for cusp-caustic violations caused by galactic-sc ale struc- 
tures. Applying the same statistical argument as 1x7 et alJ 
iooi), we conclude that the chance of observing such a vi- 
olation rate (3/5) is - 5 - 13% for a total P^°(i?cusp'^) of 
^ 20 — 30%. This can be compared with the chance of < 1% 
that we found when considering only intrinsic substructures 
(Xu et al. 2009, 2010). However, we caution that such a 
simple statistic argument does not take into account the 
full two-dimensional probability distribution in the i?cusp- 
A^ plane. 

The existing observational sample is clearly too small 
for us to reach a definitive conclusion regarding the appro- 
priateness of the ACDM model. Our main result, however, is 
that, depending on the density profiles of CDM haloes (and 
subhaloes), the line-of-sight projection effect on the flux- 
ratio anomalies of quasar images can be comparable to or 
even larger than that from intrinsic subhalos. The resulting 
cusp-violation probability from the combined effect allevi- 
ates the discrepancy between the CDM model and current 
data. New multiply-lensed four-image systems discovered in 
upcoming lensing surveys will make it possible to use the 
statistics of flux-ratio anomalies to constrain the properties 
of dark matter halo as well as the cosmogonic model. 

We end by noting that a warm dark matter cosmogony 
has a different power spectrum of density perturbations, as 
well as different density profiles for small halos compared 
to the standard CDM cosmogony. This will result in dif- 
ferent (presumably lowe r) cusp- violation probabilities (e.g., 
'Mira nda Maccial2007l ) . This possibility is worth exploring 
further in future. 
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